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We suggest that the transition to superconductivity in a single crystal of Bao.6Ko.4Bi03 with aTc — 32K, and having critical 
fields with anomalous temperature dependencies and vanishing discontinuities in specific heat and magnetic susceptibility, may 
well be an example of a fourth order (in Ehrenfest's sense) phase transition. We have derived a free energy functional for a 
fourth order transition and calculated (for the temperature range Tc/2 < T ~ Tc) the temperature dependence of the critical 
fields. We find Hci{T) cx (1 - T/Tcf, Ho{T) oc (1 - T/Tcf and Hc2iT) cx (1 - T/Tc^ in general agreement with experiments. 
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In the Ehrenfest classification of phase transitions [Q , 
an n*'' order transition is described by continuous deriva- 
tives, with respect to temperature and a mechanical vari- 
able (for example a magnetic field or pressure) up to or- 
der (n — 1). The n*^ order derivatives are discontinuous. 
So far however, only first and second order transitions 
have been observed. There are no known examples of 
transitions higher in order than two. 

We report below, what appears to be an example 
1^ of a fourth order phase transition. In the course 
of measuring ||^ the magnetization of superconducting 
Bag. eKo. 46103 (BKBO) Q , as a function of magnetic 
field (up to 27T) and temperature (1.3K to 350K) we 
were surprised to find no evidence of a discontinuity in 
the magnetic susceptibility. While this was an anomalous 
property, it was congruent with the other mystery about 
BKBO, that there is no discontinuity in specific heat ei- 
ther @ at Tc- Since in a second order phase transition, 
the boundary between the normal and superconducting 
phases satisfies, 



dH, 



c2 



dT 



AC 



(1) 



with both AC and Ax vanishing, a question arises con- 
cerning the order of this transition. 

The answer is provided by the thermodynamic crit- 
ical field Hq{T). Since in the superconducting state 
1^, the thermodynamic critical field, Hq{T), is given by 
{0 < H < Hc2), / M • dH = -H§/8Tr. This is the free 
energy of the superconducting state that is derived from 
the experimentally determined M(II,T); thus, in case of a 
second order phase transition, should have the tempera- 
ture dependence F{T) = -H^/8tt cx -{1-T/Tc)^. That 
is, Ho{T) would be linear in (1— T/Tc), apart from critical 
fluctuation effects which lead to a divergent specific heat. 
As shown in Fig. |, with = 32 K HoiT) oc (l-T/Tj^. 
Since for an n*^ order phase transition, the critical field 
has an exponent of n/2, the transition here must be of 
fourth order in the sense of Ehrenfest. 

Further support for this assertion comes from the tem- 
perature dependence of other critical fields. In particu- 
lar, we find experimentally that the lower critical field, 
the field which separates the Meissner state (no flux in 
the sample) from the Abrikosov state (partial flux pen- 
etration in the form of a vortex lattice), depends ||] 
on temperature as Hci{T) cx (1 — T/Tc)^ as shown in 
Fig. ||. The upper critical field, which separates the 
Abrikosov state from the normal state is measured to be 
Hc2iT) oc (1 - T/Tc)^-^ as shown in Fig. | This fact 
leads to an anomalous result, specific to this higher order 
phase transition. For a BCS superconductor, the ratio 
= Hc2{T) / Hci{T) is a constant. Here it diverges ap- 
proximately as (1 — T/Tc)^"^ ■ Both of the critical fields 
are inversely proportional to squares of the two length 
scales in the problem, the London penetration length A, 



which controls the flux penetration and therefore the size 
of a vortex, and the superconducting coherence length 
^ which determines the stiffness of the local density of 
the superconducting electrons. In a BCS superconduc- 
tor, these length scales are identical in their temperature 
dependence. To our knowledge there is no fundamental 
reason why k should be a constant. 

In the following we derive a free energy functional 
which describes the properties of a fourth order phase 
transition. Once we include the interaction of the su- 
perconductor with the magnetic field in the usual gauge 
invariant form, we also can derive the temperature de- 
pendencies of the critical fields. These results are in full 
accord with the experiments. 

The free energy is derived following the requirement of 
a fourth order phase transition, viz. F(T) = — /o(l — 
T/Tc)* as a function of temperature. The free energy is in 
the spirit of a Ginzburg-Landau (GL) functional which is 
minimized with respect to the complex order parameter 
ijj = Ae"^. The value at the minimum then is the ther- 
modynamic free energy. The first two terms are self evi- 
dent. Indeed it is important that the terms proportional 
to \tp\'^ and \tp\'^ be not present. The form of the spatial 
gradient term is also determined by the same considera- 
tions. The term below is the one with the lowest power 
of gradients. Higher power of gradients such as jV'V^V'l^ 
and jV'^V'P are possible but they contribute higher order 
nonlinear contributions of the magnetic field and there- 
fore are unnecessary for a stability analysis. They can be 
included for effects nonlinear in the magnetic field. The 
free energy functional appears as. 



Fiv(^, T) = a\i^f + 6|V|' + ciV'VV'l' 



(2) 



Here a ~ ao{T /T^ — 1) and b and c are positive constants. 

The minimum of this free energy corresponds to an 
order parameter amplitude A(T) oc (1 - T/TcY/"^. The 
specific heat is expected to be Civ(T) oc (1 — T/Tc)^, and 
X = oc (1 — H/Hc2)^. The thermodynamic discon- 
tinuities are in the fourth derivative of the free energy 
(or in the second derivative of specific heat as a func- 
tion of temperature or the second derivative of the mag- 
netic susceptibility x with respect to the magnetic field). 
We see that in the common thermodynamic observables, 
there are no discontinuities, as seen in the experiments. 
It is conceivable that broad transitions that have been 
observed in the past, instead of being recognized as can- 
didates for a higher order phase transition were forcibly 
squeezed into a second order framework. The Ehrenfest 
relation appropriate for a IV order phase transition is 



dH 
'dT 



A^ 



(3) 



In the presence of a magnetic field the gradient term 
transforms as V (V -I- ^A), where A is the vector 
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potential. Here (j)o is the flux quantum, <j)o = h/2e = 
2 X IQ^^^T -m? . Thus Eq. (2), as always, is the basis for a 
study of both spatial thermodynamic fluctuations as well 
as magnetic field effects. We note that the penetration 
depth for a magnetic field, the coefficient of the in the 
generalized Eq. (2), diverges as 

\-^{T) = ^^^cA^ oc (1 - T/T,f (4) 
K 

This too is in agreement with experiments, not as a di- 
rect measurement but that of iJci(T) oc (po/X^. This is 
shown in Fig. ^ with Hci{T) plotted as a function of 
temperature. Here the data on the lower critical field 
is limited by its size at temperatures close to T^. At 
low temperatures (less than Tc/2), Hci{T) behaves lin- 
ear in T and has the right intercept at Tc. The spatial 
fluctuations of the order parameter are still governed by 
^2 = c/a oc (1 - T/Tc)-^. For Hc2{T) we recognize that 
4>o/^^ — Hc2{T). Experimentally, as shown in Fig. ||, 
the exponent is nearly one. 

The proposal here rests on several critical assump- 
tions. For example, Graebner et al |^ have reported 
a very small specific heat discontinuity. The reported 
discontinuity is in fact anomalously small and roughly 
of the size of their experimental uncertainly. To esti- 
mate the expected H discontinuity, consider the specific 
heat results in ref B . The high temperature limit of 
the specific heat can be described by C(T) — + [iT"^ 
with 7 ~ 150mJ /moleK'^ as the electronic contribution 
to C{T). This large 7 puts BKBO in the category of 
heavy fermion compounds and the expected AC (of the 
order of ^Tc should be nearly 5J/moleK, considerably 
more (by a factor of 10^) than the experimental uncer- 
tainty and the reported value in ref Moreover, Hund- 
ley et al 1^ find that at low temperatures (T < Tc/S), the 
linear term in C{T) disappears. The specific heat then is 
given hy C = (3 T^, where > fi. But this larger (i may 
well be due to the presence of nodes in the putative en- 
ergy gap at the Fermi surface. For example point nodes 
in the energy gap give rise to a C oc T'^ augmenting the 
well known phonon contribution with the same power. 

Another basis for the suggestion here is the tempera- 
ture dependence of the lower critical field Hci{T). The 
cubic temperature dependence here is in contrast to the 
results of Grader et al ||] where Hci{T) is linear, as ex- 
pected for a second order BCS superconductor. However, 
closer inspection reveals that the Hd (T) values of Hall et 
al 1^ (the values used in this analysis) are at low temper- 
atures T < Tc/2 in agreement with the results of Grader 
et al who employed in their study high quality mi- 
crocrystals to eliminate spurious effects associated with 
sample inhomogeneities. The values given in ref for 
T > Tc/2, while in general agreement with ref can be 
seen to follow a straight line but extrapolate to a smaller 
Tc ~ 27K. When the zero field ^ 32 K is included. 



it is impossible to avoid a curvature in the temperature 
dependence of Hci{T). 

We note, in passing, that the relation Hq^ ~ HciHc2 is 
still valid. The consequences, near of a divergent A are 
more curious. For example the central result that the flux 
expulsion happens more slowly in the mixed state is clear. 
That the vortex lattice appears more slowly and therefore 
the irreversibility field is smaller is less obvious. Other 
questions such as the symmetry of the vortex lattice are 
currently under study and will be reported later. Sim- 
ilarly, the surface energy of a normal-superconducting 
(N-S) domain wall, is negative and proportional to A and 
therefore larger than in a BCS case. It may well engender 
a more inhomogeneous ground state at Hc2- A numeri- 
cal analysis of these questions is in progress and will be 
reported later. 

It is important to note that the thermodynamic be- 
havior changes for T < Tc/2. This is clearly seen in 
several independent measurements, for example specific 
heat and critical fields. The discussion here is focussed 
on the order of the transition from the normal state and 
is therefore limited to the vicinity of Tc. However, a mi- 
croscopic theory which might attempt to derive Eq. (2) 
will also have to include an explanation of this crossover 
behavior and possible existence of point nodes in the en- 
ergy gap. 

It might also contain an explanation of why the free 
energy does not contain terms such as and '5''. At 
present we can only speculate about a microscopic theory. 
In a sense, this question is equivalent to the seemingly 
deeper question: Why is the transition of order IV? In 
this paper, we have focussed on the properties of a IV 
order phase transition, but let's speculate: for instance, 
the BCS/GL theory contains an overall factor of density 
of states at the Fermi surface. Suppose, as discussed in 
ref. H, the density of states A^(0) = for T > and 
iV(0) oc for T < Tc. This would be a transition 
from an insulator to a superconductor, the free energy 
for p = 1 would not have a ^E"^ term and the order of the 
transition would be HI. For p=2, the transition would be 
of order IV. 

Now, addressing the relationship between fluctuations 
as developed for a II order phase transition and the 
framework: in a second order phase transition, including 
the critical effects, one might view the free energy as de- 
pending on temperature as Fo{T) = -/o(l - r/rc)^^-") 
Thus the small quantity a is calculated by pseudo- 
perturbative schemes (such as Gaussian approximation 
or some version of renormalization group). It is clear, 
however, that a value of a = —1 or — 2 is essentially 
beyond the realm of a perturbative approach. If the free 
energy exponent is significantly different from 2, then the 
unperturbed ground state could be a transition of order 
corresponding to the nearest integer, about which a cal- 
culation of fluctuations could be done in the future. 

The conclusions presented here are the first part of a 
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work in progress. Wc arc currently working on deter- 
mining (1) the magnitude of the fluctuations, and (2) 
whether there is an upper critical dimension and, if so, 
what it is. These, and other points of interest, will be 
presented in forthcoming publications. 

In summary then, we have analyzed the thermody- 
namic properties of the superconducting phase transition 
in Bao.eKo. 46103 (BKBO). The absence of a discontinu- 
ity in specific heat and magnetic susceptibility, on trans- 
forming from the normal to superconducting state, shows 
that the phase transition cannot be of second order. The 
temperature dependence of the thermodynamic critical 
field shows that the transition is fourth order. The con- 
clusions about other critical fields, derived from a free 
energy developed for a fourth order phase transition, i.e. 
Hci{T) oc (1 - T/Tcf and Hc2(T) oc (1 - T/Tc) are in 
accord with the experiments. 
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which is derived for a BCS transition, Kwok et al estimate 
ACp = 120mJ/moleK, a much smaller value than our 
estimate, but a value that is still sufficiently large to be 
clearly visible in their Cp measurements. If our proposal 
here is correct then the basis of this estimate is question- 
able, because our k is temperature dependent. 
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